Isotropic three-dimensional left-handed meta-materials 
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We investigate three-dimensional left-handed and related meta-materials based on a fully symmetric multi- 
gap single-ring SRR design and crossing continuous wires. We demonstrate isotropic transmission properties 
of a SRR-only meta-material and the corresponding left-handed material which possesses a negative effective 
index of refraction due to simultaneously negative effective permeability and permittivity. Minor deviations 
from complete isotropy are due to the finite thickness of the meta-material. 

PACS numbers: 41.20Jb, 42.25.Bs, 42.70.Qs, 73.20.Mf 



The realization of a perfect lensA and other applica- 
tions of negative refraction requires the fabrication of three- 
dimensional, homogeneous, isotropic left-handed materials 2 
(LHM) with simultaneously negative permittivity e and mag- 
netic permeability fi. So far, no such materials exist, nei- 
ther in nature nor in the laboratory. Today's available LHM 
structures, based on the periodic arrangement of split ring 
resonators 3 (SRR) and continuous metallic wires 4 , are only 
one dimensional 5 - 67 (Id), supporting left-handed properties 
only for propagation with fixed polarization in one direction, 
or two dimensional 8 - 910 (2d), where propagation in two di- 
rections with fixed polarization or one direction with arbitrary 
polarization is possible. Earlier attempts to designs at least 
isotropic SRR 11 where lacking the symmetry of SRR and unit 
cell and required individual tuning of the parameters in the 
different spatial directions. 

In this paper, we propose a three dimensional (3d) isotropic 
LHM design which allows left-handed behavior for any direc- 
tion of propagation and any polarization of the electromag- 
netic wave. Using numerical transfer matrix simulations, we 
verify the isotropic transmission properties of the proposed 
structures. Our data shows excellent agreement with results 
expected for a homogeneous slab with the corresponding neg- 
ative e and /i. 

Our meta-materials are defined as a 3d periodic continua- 
tion of a single rectangular unit cell, consisting of SRRs and 
continuous wires. The sample is a slab of meta-material with 
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FIG. 1: (Color online) The design of a fully symmetric unit cell for 
an one unit cell thick slab of an isotropic SRR (a) and a left-handed 
(b) meta-material. The interfaces are parallel to the left and right 
SRR. The metal of the 4-gap SRR (red/dark gray) and the continuous 
wires (green/medium gray) is silver using a Drude-model permittiv- 
ity around 1 THz. The SRR-gaps are filled with a high-constant di- 
electric (light gray) with a relative permittivity e gap = 300 to lower 
the magnetic resonance frequency. 



FIG. 2: (Color online) The retrieved real (solid line) and imaginary 
(dashed line) part of the effective e c s (u) and /i c ff (ui) of the homo- 
geneous medium approximation for normal incidence ($ = 0) to a 
single unit cell layer of the SRR and LHM meta-material as a func- 
tion of the dimensionless frequency u — 2nLf/c. 



a finite thickness of an integral number of unit cells and infi- 
nite extent in the perpendicular directions. The two surfaces 
of the slab are parallel to any face of the unit cell. An inci- 
dent electromagnetic plane wave with wave vector k can be 
characterized by two angles: the incidence angle {) £ [0, 7r/2) 
between k and the surface normal n of the sample, and the 
angle <f> £ (— tt, it] between the projection of k into and some 
chosen edge of the unit cell inside the surface plane of the 
sample. The frequency of the incident wave is chosen such 
that the vacuum wavelength is approximately 10 times larger 
than the linear size of the unit cell and we expect effective 
medium behavior. 

To achieve an isotropic meta-material it is wise to start with 
a cubic unit cell. The SRR and the continuous wire provide a 
resonant negative // and a negative plasmonic e, respectively. 
This occurs only when there is a component of the magnetic 
field normal to the SRR plane and a component of the elec- 
tric field parallel to the wire. Therefore we have to use one 
pair of SRR and wire for each of the three spatial directions. 
The idea is that if the interaction between the SRR and wire in 
one direction, and the mutual interaction in between the con- 
stituents in the different directions is negligible, we can, in 
first order, control the three diagonal elements of the permit- 
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tivity and permeability tensors independently. Previous results 
and extensive numerical simulations have further provided a 
set of design criteria to attain the best possible isotropy: 

(i) To avoid the coupling of the electric field to the magnetic 



resonance 



12.13 



of the circular current in the SRR we have to 
provide mirror symmetry of the SRR plain with respect to the 
direction of the electric field. Otherwise the electric resonant 
response of the magnetic resonance or the emerging cross- 
polarization terms in more-dimensional media may destroy 
the desired left-handed behavior. The problem can be avoided 
in Id and certain 2d (Id propagation with arbitrary polariza- 
tion, 2d propagation with fixed polarization) meta-materials. 
However, for an isotropic material with arbitrary direction and 
polarization of the incident wave this can only be guaranteed 
by an inversion symmetric SRR design. 

(ii) From an isotropic, homogeneous effective medium we do 
not expect any cross-polarization scattering amplitudes. Sup- 
pression of the cross-polarization terms requires, besides an 
isotropic material discretization 14 , that the inversion symme- 
try of the unit cell is not simultaneously broken in both direc- 
tions perpendicular to the direction of propagation. Therefore 
we have to center the SRRs on the faces of the unit cell and 
the wires with respect to the positions of the SRRs. 

(Hi) For the finite slab a description in terms of an homoge- 
neous effective medium is only possible if the isotropy, ie. the 
inversion symmetry of the slab, in the direction of the finite 
dimension is preservedi^A&. As a consequence, we have to 
terminate the second surface of the slab by a repetition of the 
opposite surface. This is illustrated in Fig.[^for a sample of 
one unit cell thickness. 

(iv) To avoid effects of the periodicity it is imperative to have 
the magnetic resonance frequency well below the first period- 
icity band-gap. But even in this case we observe artifacts^ in 
the effective medium response, like resonance/anti-resonance 
coupling and negative imaginary parts. These issues have 
been addressed in detail in Ref. 15. Further we have to make 
sure that the magnetic resonance frequency is below the ef- 
fective plasma frequency in the LHM, which is reduced com- 
pared to the plasma frequency of the wires by the additional 
electric (cut-wire) response of the SRR'&ii 

(v) In the LHM, the best position for the wires has been found 
to be aligned with the middle of the SRR in the center between 
a SRR and its periodic continuation. This design minimizes 
the disturbance between SRRs and wires. Additionally, it is 
favorable to put the edges of the SRRs not too close to one 
another to avoid capacitive coupling across adjacent gaps. 

To meet the symmetry requirements of the SRR we deploy 
a symmetric single-ring design with the gap distributed over 
all four sides as shown if Fig.^ The 4-gap single-ring SRR 
has been proven as a simpler, symmetric and working alterna- 
tive to the conventional design using nested rings with a single 
gap each 20 . The gap-width in the 4-gap SRR has to be signif- 
icantly reduced in comparison to a 1-gap design since due to 
the receded capacitance of the serial capacitors the resonance 
frequency of the SRR would increase. Alternatively, we could 
increase the dielectric constant inside the gaps to move the 




).4 0.45 0.5 0.55 0.6 0.65 



4 45 0.5 0.55 0.6 0.65 0.7 



FIG. 3: (Color online) Transmission spectra T — \t(uj)\ 2 for a single 
layer of unit cells of the SRR and the LHM meta-material for various 
angles of incidence (#) and polarizations (<tj, TE and TM mode). 
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FIG. 4: (Color online) Analytic transmission spectra T = |t(u;)| 2 of 
the TE and TM mode for a single layer of homogeneous unit cells of 
the SRR and LHM material for two angles of incidence. For fi(u) 
and s(lo) we used the retrieved values for normal incidence (Fig. [2}- 



resonance frequency down. In the simulations presented in 
this paper we chose the latter strategy. 

For the numerical simulations we employ a lattice transfer- 
matrix technique as described in^iS to calculate the complex 
scattering amplitudes of the sample. The finite discretization 
in conjunction with periodic boundary conditions does only 
allow a finite set of values for the component of the wavevec- 
tor parallel to the surface, hence only a finite set of inci- 
dence directions. As a consequence, to calculated the scatter- 
ing amplitudes for an incident plane wave under an arbitrary 
incidence angle $ we have to impose more general, quasi- 
periodic boundary conditions to the unit cell in both direc- 
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FIG. 5: (Color online) Transmission spectra T = \t(uj)\ 2 for a 10 
unit cell thick slab of the SRR and the LHM meta-material for vari- 
ous angles of incidence ($) and polarizations (<f>, TE and TM mode). 
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FIG. 6: (Color online) Same as Fig.|4]but for 10 unit cell thickness. 



tions perpendicular to the stratification direction. This means 
the electromagnetic fields at either side of the unit cell are 
identified up to a phase e l ® mL where m enumerates the per- 
pendicular directions, Q m S [0, 2ir/L) are the Bloch mo- 
menta parallel to the surface, and L is the size of the unit 
cell. Since the parallel momentum is preserved across the sur- 
faces, we can calculate the Q m from the angles of incidence, 
Qi = (fc|| mod27r/L) cos0, Qi = (fc|| mod27r/i) sin</>. 

We calculate the transmission and reflection amplitudes for 
a slab of the isotropic SRR and the corresponding isotropic 
LHM meta-material of one and ten unit cells thickness. The 
metals are implemented using a Drude-model permittivity for 
silver around 1 THz. The unit cell size is chosen accordingly 
to facilitate a resonance frequency in the same region. All 
results are given in terms of the dimensionless frequency u> — 



2irfL/c since up to the low terahertz region the scaling of the 
meta-materials resonance frequencies is linear and virtually 
independent on the permittivity of the metals. 

Figure|2]shows the retrieved e(uj) and fi(uj) of the effective 
homogeneous medium approximation of the SRR and LHM 
slab of one unit cell thickness for normal incidence, proving 
the presence of the magnetic resonance featuring /i < in 
the SRR and the left-handed bandpass with e < 0, fi < 
in the LHM. The retrieved e(u>) and are, in good ap- 
proximation, independent on the thickness of the slab. The 
clearly visible resonance/anti-resonance coupling and the oc- 
currence of negative imaginary parts are artifacts arising from 
the periodicity 15 which could be reduced by moving the mag- 
netic resonance to sufficiently low frequency. 

Figures |3 and |3 show the simulated transmission |£(w)| 2 
around the magnetic resonance for the one and ten unit cell 
thick slab, respectively, for different incident angles $ and 
polarizations. Because of the unit cell's symmetry it is suf- 
ficient to consider </> G [0,7r/4] for the TE and TM mode. 
Despite the square shape of the SRR, we find absolute inde- 
pendence of the scattering amplitudes (reflection not shown) 
on the orientation <f> of the incidence plane for both, TE and 
TM mode, arbitrary d, SRR and LHM meta-material slabs 
of any thickness. This finding is important for experimen- 
tal realizations as it demonstrates that a circular shape of the 
SRRs is not required for isotropic results. As expected for 
an homogeneous slab, the transmission spectra depend on the 
incidence angle -d, and, for non-normal incidence, also dif- 
fer for the TE and TM polarization. In accordance with the 
retrieved e and /i, the SRR spectra show a transmission dip 
at the magnetic resonance frequency « 0.54, which widens 
into a stop-band for the longer slab, while the LHM spectra 
show the corresponding transmission peak, developing into 
a well-defined LH passband for the longer system. This be- 
havior prevails from normal to nearly parallel incidence. Sur- 
prisingly, for off-normal incidence a second spectral feature 
emerges for both, SRR and LHM meta-material at around 
u) w 0.61 but only in the TE polarization, which gets less 
accentuated for the longer slab. However, this turns out to 
be an expected feature for an isotropic, homogeneous slab 
and has a simple explanation: The analytic transmission am- 
plitude has the form t^ 1 = cos qd — (i/2)(£ + 1/0 singe?, 
where q is the component of the wavevector perpendicular to 
the surface inside the slab, k the corresponding component 
in vacuum, and ( = fJ,k/q for the TE and ( = ek/q for the 
TM mode, what becomes the impedance of the slab for the 
case of normal incidence, q is given by the dispersion relation 
q 2 + k 2 — fj,eoj 2 /c 2 — 0. Whenever /j,(uj) for the TE mode 
or e(uj) for the TM mode becomes zero, the prefactor of the 
sine in t^ 1 diverges. For normal incidence, ie. k\\ = 0, this di- 
vergence is compensated by the vanishing q oc (/ie) 1 / 2 — » 0. 
For z9 > 0, however, fcy ^ renders the sine non-zero such 
that the divergence of the prefactor leads to a dip in the trans- 
mission t. The observed dip for the TE mode in the simula- 
tions indeed clearly corresponds to the zero of the permeabil- 
ity present in Fig.|2] further supporting the isotropic behavior 
of the meta-material slabs. 

Figures [4] and [6] show the transmission spectra analyti- 
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cally calculated for homogeneous slabs of the corresponding 
lengths, assuming homogeneous e(uj), fi(cj) taken from the 
results shown in Fig. |2] The similarity to the simulation re- 
sults, particularly for the SRR case, is striking, again confirm- 
ing the isotropic behavior of the meta-material slabs. Slight 
deviation of the simulated transmission for both, SRR and 
LHM, can be attributed to two principal deficiencies of our 
estimation of e(w) and First, the artifacts arising from 

the periodicity will introduce an explicit dependence of the in- 
cident angle $ since the unit cell is rectangular (and we expect 
the band structure to look different for various directions). 
Second, due to the symmetric surface termination the relevant 
periodicity for the retrieval shown in Fig.|2]is one mesh step 
longer than the period in the perpendicular directions. These 
effects will become negligible if we move the magnetic reso- 
nance to sufficiently low frequency while maintaining the size 
of the unit cell. 

The LHM transmission spectra reveal an additional, more 
serious problem. The finite thickness of the slab leads to an 
essential asymmetry in the electric response between the di- 
rections parallel and normal to the surface. The electric field 
in the normal direction sees a finite, instead of a continuous 
wire, which changes the electric response, hence e(uo) in that 
direction, leading to an deformed spectrum for large -d. This 
problem only affects the TM mode for oblique incidence. As 
is apparent from Fig. [5] within the LH passband we can expect 
this issue to become less important with increasing thickness 
of the slab. However, in the transmission gaps the near unity 
reflection is dominated by the few surface layers of the slab, 
thus experiencing the termination of the finite wires for any 
system length. 

We are aware that the geometry of the SRR with dielectric- 
filled gaps shown in Figgis not very convenient for experi- 



mental realizations. However, the issues discussed in this pa- 
per are generic and do not require this particular realization of 
the SRRs. Experimentally more realistic SRRs that can oper- 
ate at very low frequency not sacrificing the symmetry can be 
designed based on multi-layer 4-gap structures or two-layer 
opposing conventional SRRs 13 . For low frequencies also the 
wires may be thinner, such that the crossing can easily re- 
placed by centered wires which just touch each other in the 
center of the unit cell. 

In conclusion, we presented a design for an LHM and 
SRR-only meta-material based on a fully symmetric unit 
cell. Transfermatrix simulations of finite thickness slabs for 
oblique incident plain waves revealed almost isotropic trans- 
mission spectra. Further, the simultaneously negative e and 
/i has been confirmed for the LHM slab. We emphasized 
the importance of symmetry issues in the design of isotropic 
meta-materials. Especially the asymmetry introduced by the 
gaps in the conventional SRR design has to be avoided. Addi- 
tionally, the magnetic resonance frequency has to be chosen 
much lower than the first band gap arising from the meta- 
materials intrinsic periodicity to avoid interference with the 
band structure and the breakdown of the effective medium be- 
havior. Minor deviations from the isotropic behavior which 
occur particularly in the LHM slab have been attributed to the 
finite thickness of the meta-material which modifies the elec- 
tric response in the short direction. Most of the imperfections 
become unimportant for sufficiently long systems. 
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